Spin-orbit gap of graphene: First-principles calculations 
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Even though graphene is a low energy system consisting of the two dimensional honeycomb lattice 
of carbon atoms, its quasi-particle excitations are fully described by the 2+1 dimensional relativistic 
Dirac equation. In this paper we show that while the spin-orbit interaction in graphene is of the 
order of imeV, it opens up a gap of the order of W'^meV at the Dirac points. We present the first 
principle calculation of the spin-orbit gap, and explain the behavior in terms of a simple tight-binding 
model. Our result also shows that the recently predicted quantum spin Hall effect in graphene can 
only occur at unrealistically low temperature. 
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PACS numbers: 73.43.-f, 71.70.Ej, 73.21.-b 

Recently, electronic properties of graphene, a single- 
layer graphite sheet, attracted great interests both the- 
oretically and experimentally. The key difference of 
graphene compared with most other two-dimensional ma- 
terials is the linear energy spectrum around two nodal 
points in Brillouin zone, which makes the low energy dy- 
namics of electrons in this system equivalent to that of 
relativistic fermions, as described by the massless Dirac 
equation [l[ . The two sublattices in graphene honeycomb 
lattice playthe role of pseudo-spin degrees of freedom. 
In Ref. 0,13, the quantum Hall effect in graphene is ob- 
served, which shows non-conventional quantization rule 
an = ^{'^n + G Z. Such an "abnormal" quan- 
tum Hall effect agrees with the theoretical calculation 
based on massless Dirac equation under external mag- 
netic field [1, @, @|j and can be considered as a conse- 
quence of the chiral anomaly in two dimensional massless 
fermions. Moreover, recent experiment on low field mag- 
netoresistance Q shows that graphene remains metallic 
under temperature as low as T = AK, which confirms 
that any possible gap opened at the Dirac cones cannot 
be larger than ksT ~ 0.34meV. 

Nevertheless, it has been proposed that a small gap 
can open on the two Dirac points of graphene due to 
spin-orbital coupling (SOC) [8|, which at the same time 
makes the system a spin-Hall insulator P| with quan- 
tized spin Hall conductance. Physically, this proposal is 
a spinful version of Haldane's model for quantum Hall 
effect without magnetic field [Tc| , in which a spin depen- 
dent next-nearest-neighbor hopping term is introduced 
to induce opposite mass terms for the two Dirac cones. 
Ref. Q estimates the spin-orbit gap in graphene to be 
2AK. In this paper we provide systematic calculations 
of the spin-orbital gap in graphene by both first-principle 
calculation and tight-binding model, which shows that 
the actual gap is much smaller compared to the crude 
estimate given in Ref. it explains the (nearly) gap- 
lessness observed in experiments and also defines a much 



more narrow temperature range for the quantum spin 
Hall effect to be observed. 

The sp^ hybridization of the 2s orbital and two 2p or- 
bitals of carbon atom creates the a bonds to form the 
honeycomb lattice of graphene which is bipartite with 
two carbon atoms in one unit cell. The n band consist- 
ing of the remaining 2p orbitals controls the low energy 
physics of graphene and makes it a semi-metal. One can 
describe the tt and a electrons by two tight binding (TB) 
Hamiltonians separately, which in momentum space is a 
2x2 matrix H^(/c) for tt band, and a 6 x 6 matrix ll^(k) 
for cr band [11]. If the spin degeneracy of electrons is 
taken into account, the dimension of these two matrices 
are doubled. The diagonal entities of the matrices are the 
on-site energies of different orbitals and the off-diagonal 
entities are the possible hopping between different sub- 
lattices. 

The SOC is a relativistic effect described by the Hamil- 
tonian with the form Ha ■ {W x p)/{Am?c^) ^ L ■ a. a 
is the Pauli matrix. For a single carbon atom, there is 
no SOC between 2s and 2p orbitals due to their differ- 
ent azimuthal quantum number, and the SOC only exists 
among the 2p orbitals. Its magnitude can be estimated 
of order AmeV by directly computing the overlap inte- 
gral of SOC between 2pz and 2pa; orbitals. Note that the 
SOC changes the magnetic quantum number accompa- 
nied with the spin-flip of electrons, hence no SOC exists 
between the same atomic orbital. 

For the graphene, only the SOC in the normal direc- 
tion with the form L^cTz has nonzero contribution due to 
the reflection symmetry with respect to the lattice plane. 
Even this term vanishes for the tt orbitals between near- 
est neighbors (NNs), since there is an additional vertical 
reflection plane along the nearest neighbor bond. Under 
the mirror reflection to this plane, the 2pz wavefunctions 
of the adjacent atoms are unchanged, however the an- 
gular momentum changes its sign, hence the matrix 
element of LzCJz between NNs vanishes. This is different 
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from the carbon nanotube[l^ where the curvature effect 
could provide a SOC between 2p orbitals of NNs although 
it is still vanishing for large tube radius. Thus, to realize 
the SOC effect of tt band within the NNs approximation 
we need the aid of cr-band. This process turns out to be 
a second order one, which is three orders of magnitude 
smaller than ^o- On the other hand, the SOC can act 
directly within the cr-band, it will open a gap at some 
degenerate points with the same order as £_o ~ 'imeV. 

The SOC mixes the tt and cr bands and the total Hamil- 
tonian reads 



H = 



a. T 



(1) 



Here, H^r and Ho- should be enlarged to be 4 x 4 and 
12 X 12 matrices by the spin indices, respectively. The 
SOC term T bridging the tt and cr bands is a 4 x 12 
matrix of order f^nd its explicit form is not important 
at present, and will be given later. The wavevector k is 
omitted for simplicity hereinafter because it is always a 
good quantum number. 

Since we are concerned with the low energy physics, 
an effective 7r-band model with SOC derived from the 
original Hamiltonian Eq. ^ is more advantageous. For 
this purpose, one can perform a canonical transformation 



H - 
S = 



Hs = e-'^He" 
M 
-Mt 



where M should satisfy 



(2) 



(3) 



so that Hs is block diagonal up to order ^q. Clearly M 
is also a 4 X 12 matrix. The effective Hamiltonian i?e// 
is then extracted from the diagonal part of Hs as 



iIe//«H,-i(TMt+MTt) 



(4) 



The second term is just the effective SOC for the vr band 
electrons. 

The matrix M can be calculated iteratively through 
Eq. © 



M = TH-i+H^TH„ 



(5) 



Around the Dirac points, the spectrum of H,r is close to 
zero measured from the on-site potential of 2p orbital, 
while that of Ho- is of order several eV , hence we can 
take M w TH^^ approximately. The effective SOC of 
TT-band then reads 



(6) 



whose magnitude is roughly estimated as ~ j^oP/^ 
with A being of the order of the energy difference at the 



Dirac points between tt and cr bands, is of the order 
l{)~^meV , since A is of order eV . 

So far we have not used the explicit form of Hg. and 
T in the above discussions. To derive and the SOC 
analytically, we need more details of Ho- and T. can 
be written as 



E S 

St E 



I, 



(7) 



where I is the identity matrix for the spin degrees of 
freedom. The matrix E represents the onsite energy of 
different atomic orbitals, which can be written as 




E= I 
A^ 



(8) 



if we arrange the three sp^ hybridized orbitals in the se- 
quence of {2py, 2px, 2s}. Here, A^ is the energy difference 
£2s — £2p between the 2s and 2p orbitals. S describes the 
hopping between the two sublattices in the momentum 
space. To give its exact form, we first consider the hop- 
ping between the two adjacent atoms in real space, which 
can also be described by a 3 x 3 matrix. Suppose the two 
adjacent atoms are placed on the horizontal x-axis, i.e., 
the bond angle is zero, this hopping matrix can be writ- 
ten as following 



^Pp-K 







V^pcr ^spa 

^SpG ^SS(T 



(9) 



One can obtain the hopping matrix S(6') for arbitrary 
bond angle 6* by a rotation R(^?) in the xy plane as 



I](6i) = Rt(6')I]oR(0) 



R(0) 



cos 9 — sin 6 
sin 6 cos 6 
1 



(10) 



The parameters V^pjr, Vppa, Vspa and Vssa correspond 
the cr or TT bonds formed by 2s and 2p orbitals, whose 
empirical value can be found in textbooks, for example, 
Ref. [ll|. Note that we do not consider the wavefunction 
overlap matrix in our TBA scheme for the sake of sim- 
plicity. Then the hopping matrix in the momentum space 
reads 



(11) 



where, da with = 1, 2, 3 are the bond vectors connect- 
ing the carbon atom and its three nearest neighbors and 
9a is the angle between da and x-axis. 

For T, as we have described above the spin flip on the 
same atom only takes place between the 2pz and two in- 
plane 2px^y orbitals. A straightforward calculation leads 
to the onsite spin flip 



^o(-cr2;, cry, 0), 



(12) 
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with Ux^y the Pauh matrices. Then T can be written as 



To 
To 



(13) 



Notice that there are two To terms in the above matrix 
corresponding to different sublattices. 

Since Ho^ has a large gap near the Dirac points K and 
K* , we can expect that Hcrik + K) = Hct(-/?) + o(fc), 
which means we can substitute H~^(i?) into Eq.® as a 
good approximation. Finally we get the effect Hamilto- 
nian with SOC at the low energy scale, 



H 



H 



[K] 
eff 

[K'] 



vp{kx-iky) -^icr^ 

-£.icrz vpik^-iky) 
vpik^ + iky) £,i(7z 



.(14) 



The off-diagonal terms in the above equations come from 
the well-known form of Hjr, and is just the Fermi 
velocity of tt electrons at the Dirac points. The effective 
SOC ^1 in our TBA scheme has an explicit form 



6 « |eo|'(2A,)/(9K' ) 



(15) 



Eq. (|T5|) is the key result from our tight-binding cal- 
culation. Eq. (fT4|) leads to a spectrum E[k) = 
{vpkY + Taking the values of the correspond- 



ing parameters from Ref.[ll|, one can estimate is of 
order lQ~^meV , so is the energy gap at the Dirac 
points. 

Eqs. (fT4|) are similar to those in Ref. except that 
the SOC constant is three orders of magnitude smaller 
than their estimate. We can also consider the SOC of tt 
orbitals between next nearest neighbors (NNN) which is 
not forbidden by the symmetry. In this case the electron 
moving between NNN will be accelerated by the atoms 
other than these two NNN ones which provides the cor- 
responding SOC. This will involve three center integrals, 
i.e., two orbital centers and a potential center which are 
different with each other. Generally speaking, such in- 
tegrals are very small which leads to the SOC of order 
at most lQ~'^meV by our estimate and may actually be 
smaller. 

The argument above is supported by accurate first- 
principles calculations based on density-functional the- 
ory (DFT). The relativistic electronic structure of 
graphene was calculated self-consistently by the plane- 
wave method [l3| using the relativistic fully separable 
pseudopotential in the framework of noncoUinear mag- 
netism . The exchange-correlation potential is treated 
by the local density approximation (LDA) whose validity 
for the system considered here has been demonstrated by 
many other studies. The experimental lattice parameter 
a = 2.456A is used in the calculation. The convergence of 
calculated results with respect to the number of k points 
and the cut-off energy has been carefully checked. 




FIG. 1: (color online). Calculated relativistic band structure 
of graphene. 



Fig. 1 shows the band structure of graphene, we can 
find that gap induced by SOC for a orbit is 9.0 meV at 
r point. The figure also indicates that there is a gap 
induced by SOC for tt orbit at K point, and the mag- 
nitude of the splitting gap is 0.8 x 10~^meV, which is 
in good agreement with the estimate obtained from the 
tight-binding model discussed above. Since the number 
discussed here is so small, a few notes are necessary: 1) 
the calculations are valid within LDA; 2) the numeric ac- 
curacy of present calculations reachs 10~^ meV per atom; 

3) the convergency of gap size with respect to the number 
of k points and cut-off energy is better than 1 x 10~'*meV; 

4) the Kramer doublet degeneracy can be reproduced 
down to 10~^meV. Nevertheless, it is clear that the gap 
induced by SOC at the K point is of order 10~^meV. 
Considering that graphene may be typically deposited 
on substrates, the graphene are generally strained due 
to small lattice mismatches, thus the lattice mismatch 
strain can tune the splitting gap at K point. We have 
calculated band structure for different lattice constants 
of graphene, and have found that the splitting gap in- 
creases a little with compression while the gap decreases 
with tensile strain. 

In conclusion, we provided a careful calculation on 
spin-orbit gap of graphene, which leads to the same mass 
term for the relativistic Dirac fermions in the contin- 
uum limit [11, but with a much smaller magnitude of 
the gap 10~^meV. The physical reason of the small- 
ness of the spin-orbit gap can also be understood from 
the tight-binding model as coming from the lattice C3 
symmetry, which leads to the vanishing of the leading 
order contributions. Such a small gap is consistent with 
the experimental observation of semi-metallic behavior of 
graphene. It shows that the proposed quantum spin Hall 
effect in graphene cannot be observed until temperature 
as low as T ^ 10~^K. In addition, impurity scattering 
in a disordered sample may also destabilize the effect. 
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